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Abstract 

We study the energy loss of a heavy quark propagating through the matter produced in the 
collision of two sheets of energy [1]. Even though this matter is initially far-from-equilibrium 
we find that, when written in terms of the energy density, the equilibrium expression for heavy 
quark energy loss describes most qualitative features of our results well. At later times, once a 
plasma described by viscous hydrodynamics has formed, the equilibrium expression describes 
the heavy quark energy loss quantitatively. In addition to the drag force that makes it lose energy, 
a quark moving through the out-of-equilibrium matter feels a force perpendicular to its velocity. 



The discovery that strongly coupled quark-gluon plasma (QGP) is produced in ultrarelativis- 
tic heavy ion collisions has prompted much interest in the real-time dynamics of strongly coupled 
non-Abelian plasmas. For example, heavy quark energy loss has received substantial attention. 
If one shoots a heavy quark though a non-Abelian plasma how much energy does it lose as it 
propagates? This question has been answered [2] for equilibrium plasma in strongly coupled 
TV = 4 supersymmetric Yang-Mills (SYM) theory in the large number of colors Nc limit, where 
holography permits a semiclassical description of energy loss in terms of string dynamics in 
asymptotically AdSs spacetime [3, 4, 2]. One challenge (not the only one) in using these results 
to glean qualitative insights into heavy quark energy loss in heavy ion collisions is that a heavy 
quark produced at time zero must first propagate through the initially far-from-equilibrium matter 
produced in the collision before it later plows through the expanding, cooling, near-equilibrium 
strongly coupled QGP. In this contribution we sketch our work in progress toward obtaining 
guidance for how to meet this challenge. We want a toy model in which we can reliably calcu- 
late how the energy loss rate of a heavy quark moving through the far-from-equilibrium matter 
present just after a collision compares to that in strongly coupled plasma close to equilibrium. 

We study the energy loss of a heavy quark moving through the debris produced by the col- 
lision of planar sheets of energy in strongly coupled SYM theory analyzed in Ref. [1]. The 
incident sheets of energy move at the speed of light in the +z and -z directions and collide at 
z = at time t = 0. They each have a Gaussian profile in the z direction and are translation- 
ally invariant in the two directions x± = {x\,X2} orthogonal to z- Their energy density per unit 
transverse area is ir'iNl jln^), with [i an arbitrary scale with respect to which all dimensionful 
quantities in the conformal theory that we are working in can be measured. The width a of the 
Gaussian energy-density profile of each sheet is chosen to be cr = 1/(2^). 

We study heavy quark energy loss by inserting a heavy quark moving at constant velocity 
P between the colliding sheets before the collision and calculating the force needed to keep its 
velocity constant throughout the collision. The energy density of the colliding sheets as well as 
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two sample trajectories of a quark getting sandwiched between them are shown in the left panels 
of Figs. 1 and 2. Via holography, the colliding planar sheets in SYM theory map into colliding 
planar gravitational waves in asymptotically AdSs spacetime [1]. The addition of a heavy quark 
moving at constant velocity p amounts to including a classical string attached to the boundary 
of the geometry [4] and dragging the string endpoint at constant velocity pulling the string 
through the colliding gravitational wave geometry. The force needed to maintain the velocity of 
the string endpoint yields the energy and momentum loss rates of the heavy quark [2] . 

We employ inf ailing Eddington-Finkelstein coordinates to describe the asymptotically AdSs 
geometry. Translation invariance in the x_l plane allows one to write the metric Gmn in the form 

ds^ = GuNdX^dX^ = -A dv^ + I? [e^dxl + e'^^dz^) + 2Fdvdz-2du dv/u^ (1) 

where A, B, 2 and F depend on coordinates v, z and u. Here u is the AdS radial coordinate with 
the AdS boundary located at w = 0. Lines of constant v and (z,Xj_) are radially inf ailing null 
geodesies affinely parameterized by w. At w = the time coordinate v coincides with the time 
coordinate t in SYM. For initial conditions consisting of two colliding gravitational waves, the 
functions A, B, F and S were obtained numerically in Ref. [1] by solving Einstein's equations. 
We obtain the string equations of motion from the Polyakov action 

Sf = J do-dT£p = -^^ J do-dT^f^T]''^GMN^aX^^bX^ (2) 

where Tq and L are the string tension and AdS radius and where T]ab is the worldsheet metric with 
T and cr being worldsheet coordinates chosen such that 

/ -a(T,cr) -1 \ 
r]ab = \ _^ Q I (3) 

with Qf(T, cr) to be determined. With this choice of worldsheet coordinates, which are the world- 
sheet analogue of infalling Eddington-Finkelstein coordinates, lines of constant r are infalling 
null worldsheet geodesies and the string equations of motion take the particularly simple form 

d^X^ + r^s^^d^X^ = with X^ = drX^ -\a(T,o-)d^X^ , (4) 



AB 

where F^^ are Christolfel symbols. X^ is the directional derivative along outgoing null world- 
sheet geodesies. The constraint equations are (do-X)^ = and X^ = 0. The first is a temporal 
constraint: if it is satisfied at one value of r and the other equations are satisfied, it remains satis- 
fied at all times. = is a boundary constraint: if it is satisfied at one value of cr and all other 
equations are satisfied, it is satisfied at all values of cr. Residual worldsheet diff'eomorphism in- 
variance allows us to choose u = cr. Eqs. (4) now prescribe the r-evolution of X^, as follows: 
given X^ at one r, use (4) to compute X^; with w = cr the definition of X^ becomes a = -2X^; 
with a in hand, construct djX^\ finally, determine X^ at the next r. To execute this procedure 
and obtain v(t, cr), z(t, cr) and x_l(t, cr) we first need an initial profile of the string that satisfies 
{d^Xf = and we need to impose boundary conditions at w = 0, requiring that the endpoint of 
the string atu = moves at constant velocity j3 and that X^ = remains satisfied at w = 0. After 
we have solved (4) and discovered how the string is buff'eted by the colliding sheets, we compute 
the canonical momentum flux down the string atu = 0. This flux is minus the force required 
to maintain the velocity of the string endpoint and yields the four momentum lost by the quark 
per unit time 

^ = < = ^^G^^rr - air, a-) d^r)l__^ , (5) 
2 



Figure 1: Left: The path in (t,z) of a quark moving in the x-direction with /3 = 0.5 along z = 0. Energy density of 
the coUiding sheets (in units of ju^ / (In^)) is shown by the color. Right: Instantaneous momentum loss of a heavy 
quark with two different initial conditions for the string profile (solid blue and dot-dashed red curves). The results agree 
for t > 0, illustrating that our calculation of the momentum loss after the collision is insensitive to the initial string 
profile. The dashed cyan, magenta or green curves show what the momentum loss would be in an equilibrium plasma 
with temperature Tg, Tj_ or T\\ determined from the instantaneous energy density, transverse pressure or parallel pressure. 



where A = g^Nc is the 't Hooft coupHng in the gauge theory. We choose several sets of initial 
conditions for the string and demonstrate that they all yield the same dp^ldt at suitably late times. 

The instantaneous momentum loss dpjdt experienced by the heavy quark is given by the 
solid curves in Figs. 1 and 2. These curves are the central results of our calculation. Fig. 1 
illustrates that we are initializing the string profile early enough that our results for the drag force 
are sensitive to how we do this only before the collision happens, not after the sheets collide 
> 0). In Fig. 2, the quark is at z = at ^ = but moves off' the z = plane with some nonzero 
rapidity. In this case, in the frame in which we do the calculation the ffuid is not at rest at the 
location of the quark for ^ > 0. However, in Fig. 2 we plot dpjdt in the local fluid rest frame, 
meaning that at each instant in time we transform to a frame in which the fluid at the location of 
the quark is static, but of course not in equilibrium. In particular, when the quark is moving with 
nonzero rapidity the gradient of the fluid velocity at the location of the quark is nonvanishing. 
The blue curve shows the component of dpidt in the direction opposite to the velocity vector of 
the quark in the local fluid rest frame. This is a drag force, just as in the case of motion with zero 
rapidity in Fig. 1. The red curve in Fig. 2 shows the component of dpjdt perpendicular to the 
quark velocity in the local fluid rest frame. A perpendicular force like this can arise because of the 
presence of gradients in the fluid velocity. (We define the sign of the perpendicular component 
such that dp^jdt > when \dpl^/dt\, see below, is greater than \dpz/dt\.) 

We are interested in comparing our result to expectations based upon the classic result for 
the drag force required to move a heavy quark with constant velocity through the equilibrium 
plasma in strongly coupled SYM theory, namely [2] 



dp 
dt 



(6) 

eq 
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where T is the temperature of the equilibrium plasma. Out of equilibrium, the matter does not 
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Fluid rest frame 




Figure 2: Left: The path in {t, z) of a quark with nonzero rapidity, specifically yS^ = 0.2 and /3x = 0.5. Right: Instantaneous 
momentum loss in the fluid rest frame, antiparallel to the quark velocity (solid blue) and perpendicular to it (solid red). 
Green dashed curve shows the drag force in an equilibrium plasma with temperature Te . 



have a single temperature. We can nevertheless use (6) to frame expectations for dp/dt at any 
point in spacetime, as follows. We transform to the local fluid rest frame and define three "tem- 
peratures" from the stress-energy tensor by writing it as T^y = (tt^A^^/S) diag (r^, T^i^). 
(In static equilibrium, = = T\\ = T .) We can then use each of these three "temperatures" 
in the equilibrium expression (6) and see to what degree the three curves that result bracket our 
result for the energy loss of a heavy quark in far-from-equilibrium matter. We do this in Fig. 1. 
The green dashed curve in Fig. 2 shows the force (6) in an equilibrated plasma with T = T^. In 
an equilibrated plasma, the quark only feels a drag force, opposing its velocity; the component of 
dpidt perpendicular to the velocity of the quark vanishes. The red curve in Fig. 2 arises because 
the heavy quark is propagating through out-of-equilibrium matter featuring velocity gradients. 

With the exception of the (small) component of dpjdt depicted by the red curve in Fig. 2, our 
results are rather well described by the equilibrium expression (6) with T = T^. At late times, we 
see from Figs. 1 and 2 that the drag force that we have calculated agrees quantitatively with this 
equilibrium expectation. This agreement begins at the same time that viscous hydrodynamics 
starts to describe the still far-from-isotropic fluid well [1]. At early times, when the fluid is far 
from equilibrium, (6) with T = Te does a reasonable job of characterizing our results, although 
perhaps using T = Tj^ does even better. Certainly there is no sign of any significant "extra" 
energy loss arising by virtue of being far from equilibrium. The message of our calculation seems 
to be that if we want to use (6) to learn about heavy quark energy loss in heavy ion collisions, it 
is reasonable to apply it throughout the collision, even before equilibration, defining the T that 
appears in it through the energy density. The error that one would make by treating the far-from 
equilibrium energy loss in this way is likely to be much smaller than other uncertainties. 
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